DIGITAL FILTERING METHODS AND DEVICES 

The present invention relates to digital filtering methods and devices. 
It finds an advantageous application to the filtering of digital images, 
5 in particular for the compression of images in accordance with "JPEG 2000", 

currently being drafted. 

Two techniques are notably known for implementing the filtering of 

digital signals using wavelet filters: 

- filtering by convolution, and 
10 - filtering by "lifting". 

A few elementary concepts on wavelet filters, and then on lifting, are 
given below. 

It is known that the wavelet filtering of a monodimensional signal Xj 
can be defined as: 

15 - the convolution of a low-pass filter H 0 (k) with the signal Xj, the 

signal Y n resulting from the convolution being downsampled by a factor of two: 

Y 2 n=X H 0(k).X 2 n-k (1) 
k 

and 

- the convolution of a high-pass filter Hi(k) with the signal Xj, the 
20 signal Y n resulting from the convolution being downsampled by a factor of two: 

Y 2n+1 = X! H l( k )* X 2n+1-k (2) 
k 

The result of the filtering, referred to as a wavelet filtering operation 
(Ho, Hi), produces a signal Yj, in which each even-indexed sample Y 2n is a low- 
pass coefficient and each odd-indexed odd-indexed sample Y 2n +i is a high-pass 
25 coefficient. 

The low-pass filter and the high-pass filter must satisfy certain 
conditions, referred to as perfect reconstruction conditions, so that the perfect 
reconstruction of the signal X, from the signal Yj is possible. 

If the notation of the z transform is used for the finite impulse 
30 response filters, namely 



H i (z) = XH i (k).z- k , (3) 

k 

then the perfect reconstruction condition which the filters must satisfy is: 

H 0 (z).H 1 (-z- 1 ) + H 0 (-zJ.H, (z~ 1 ) = 2 . (4) 
It is known that some wavelet filters have symmetry properties. 
5 An odd-length symmetric wavelet filter H(k) or WSS (Whole-Sample 

Symmetric) filter is defined by the following equation, referred to as WSS 
symmetry, between all the coefficients of the filter: 

Vk, H(-k) = H(k) (5) 
An even-length symmetric wavelet filter or HSS (Half-Sample 
10 Symmetric) filter is defined by the following equation, referred to as HSS 
symmetry, between all the coefficients of the filter: 

Vk, H(-k) = H(k + 1). (6) 
An even-length antisymmetric wavelet filter or HSA (Half-Sample 
Antisymmetric) filter is defined by the following equation, referred to as HSA 
1 5 symmetry, between all the coefficients of the filter: 

Vk, H(-k) = -H(k + 1). (7) 
The technique of lifting, mentioned above, consists of implementing a 
filtering step by decomposing it into an equivalent sequence of elementary 
filtering operations, referred to as lifting steps. This decomposition is referred to 
20 as lifting factorisation. For more details on lifting factorisation, reference can 
usefully be made to the article by Ingrid DAUBECHIES and Wim SWELDENS 
entitled "Factoring wavelet transforms into lifting steps", in J. Fourier Anal. 
Appl., vol. 4, No. 3, pages 247 to 269, 1998. 

It should be stated in particular that, when a filtering operation is 
25 performed in the form of lifting steps, in order to transform a signal Xj into a 
signal Yi, it is generally begun by copying all the samples of the signal X\ into 
the signal Yj (i.e. Yj = X,), and then the even-indexed samples and the odd- 
indexed samples of the signal Yj are modified alternately. 
There are two types of lifting step: 
30 - high-pass lifting steps, and 

- low-pass lifting steps. 



A high-pass lifting step consists of modifying odd-indexed samples 
by adding to them a sample function of even-indexed samples which are 
weighted or filtered: 

f 



Vn,Y 2n+1 =Y 2n+1 +R 



X a k- Y 2n+2k 

^ k 



(8) 



5 The function R is, in general terms, any approximation of the variable 

x. The function R can be a rounding operator which rounds a real value x to an 
integer (such as the closest integer), or can simply be identity: R(x) = x. 

A low-pass lifting step consists of modifying even-indexed samples 
(low-pass) by adding to them a function of odd-indexed samples which are 
1 0 weighted or filtered: 

r 



Vn,Y 2n =Y 2n+ R 



S b k- Y 2n+2k+1 



(9) 



In the case where R(x) = x, at each lifting step, each low-pass 
sample Y 2n can be expressed in the form of a linear combination of samples of 
the signal X u corresponding to a low-pass filtering of the signal X h Thus each 
15 low-pass lifting step can be associated with a low-pass filter. 

It is said that each low-pass lifting step lifts a low-pass filter (one 
which is associated with the previous low-pass lifting step) into another low- 
pass filter (the one which is associated with the current low-pass lifting step). 
The term lifting accounts for the fact that the support of the resulting low-pass 
20 filter (that is to say the number of coefficients of the filter) is greater than that of 
the corresponding low-pass filter at the previous step. The low-pass filter with 
the larger support finally obtained by lifting is called the "equivalent low-pass 
filter". 

The same applies to the high-pass filters and the lifting steps, that is 
25 to say each high-pass lifting step can be seen as the lifting of a high-pass filter 
into a high-pass filter with a larger support, referred to as the "equivalent high- 
pass filter". 

The same terminology is used in the case where the approximation 
function R is not identity. 
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Thus the lifting-based implementation of a wavelet filtering operation 
(H 0) Hi) can be seen as a sequence of operations during which, alternately and 
iteratively, low-pass filters are lifted into equivalent low-pass filters with a larger 
support, and high-pass filters are lifted into equivalent high-pass filters with a 
5 larger support, until the low-pass filter is equal to H 0 and the high-pass filter is 
equal to H<|. 

Several types of symmetry are also encountered in the lifting steps: 

- there are two types of symmetric lifting steps: 

• lifting steps with WSS symmetry, and 
10 • lifting steps with HSS symmetry; and 

- there are two types of antisymmetric lifting steps: 

• lifting steps with WSA antisymmetry, and 

• lifting steps with HSA antisymmetry. 

An odd-length symmetric lifting or WSS lifting ("Whole-Sample 
15 Symmetric lifting") step is defined by a WSS symmetry relationship between the 
filtering coefficients of the lifting step: 

a_ k =a k (10) 

An even-length symmetric lifting or HSS lifting ("Half-Sample 
Symmetric lifting") step is defined by an HSS symmetry relationship between 
20 the filtering coefficients of the lifting step: 

a„ k = a k _>, (11) 

An odd-length antisymmetric lifting or WSA lifting ("Whole-Sample 
Antisymmetric lifting") step is defined by a WSA symmetry relationship between 
the filtering coefficients of the lifting step: 
25 a 0 =0 and a_ k = -a k (12) 

An even-length antisymmetric lifting or HSA lifting ("Half-Sample 
Antisymmetric lifting") step is defined by an HSA symmetry relationship between 
the filtering coefficients of the lifting step: 

a_ k =-a k _ 1 (13) 
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A single-coefficient lifting or SC lifting ("Single-Coefficient lifting") 
step is defined by the fact that all the filtering coefficients of the lifting step 
except one (denoted a k .) are zero. 

For the low-pass lifting steps, 
5 - if k'= 0 , the lifting step is called a Right-side Single Coefficient low- 

pass lifting or RSC low-pass lifting step, and 

- if k'= -1, the lifting step is called a Left-side Single Coefficient low- 
pass lifting or LSC low-pass lifting step. 

In a similar manner, for the high-pass lifting steps, 
10 - if k'=1, the lifting step is called a Right-side Single Coefficient 

high-pass lifting or RSC high-pass lifting step, and 

- if k'=0, the lifting step is called a Left-side Single Coefficient high- 
pass lifting or LSC high-pass lifting step. 

Given that, at the present time, more and more signal compression 

15 algorithms use wavelet filtering operations for the decorrelation of the signals, it 
is very interesting to have available an effective implementation of these filtering 
operations. Various image compression algorithms use bidimensional wavelet 
transformations, which consist of the successive application of 
monodimensional transformations of the type described in the introduction. 

20 An effective implementation of wavelet-based monodimensional 

filtering is lifting-based filtering. 

It is known in fact that lifting-based filtering reduces the number of 
filtering operations (multiplications and additions) and also affords a reversible 
implementation of wavelet-based filtering, that is to say with perfect 

25 reconstruction of the signal, without loss of information, in the case of signals 
whose sample values are integer. 

In addition, the step which is the reverse of a lifting step is obtained 
in a trivial fashion by changing the sign before the approximation function R 
given in equations (8) and (9) above. This makes it possible notably to use the 

30 same circuitry in a hardware implantation corresponding to both forward and 
reverse wavelet transformations. 
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For a given pair of wavelet filters (Hq.H-i) the problem consists of 
finding a sequence of lifting steps which will implement the wavelet 
transformation. This is the problem of lifting factorisation. 

Two families of wavelet filters are known which have good signal 
5 decorrelation properties: the family WSS/WSS, abbreviated to WSS, and the 
family HSS/HSA, the notation A/B meaning that the low-pass filter H 0 has a 
symmetry A and the high-pass filter Hi has a symmetry B. 

The problem consisting of implementing WSS wavelet filters by using 
HSS lifting steps has been resolved and is described in the article by Ingrid 
10 DAUBECHIES and Wim SWELDENS cited above. 

The problem consisting of implementing HSS/HSA wavelet filters has 
been partly resolved, by using WSA lifting steps. 

However, this solution is only partial since this implementation of the 
prior art applies only to a limited set of HSS/HSA wavelet filters and not to all of 
15 these filters. 

Moreover, the publication by E. MAJANI entitled "Low-complexity 
wavelet filter design for image compression", TDA Progress Report 42-119, 
NASA Jet Propulsion Laboratory, Pasadena, CA, 15 November 1994, indicates 
that it is possible to use WSA lifting steps for lifting a pair of filters (Hq.Hj), 
20 where the two filters have the same support, into a pair of filters with a larger 
support (H , 0 ,H , 1 ). However, no indication is found therein about the way of 
implementing, by lifting steps, pairs of filters (Hq,^) where the two filters have 
the same support and where this support is greater than two. 

The purpose of the present invention is to remedy this drawback. 
25 It proposes an implementation by lifting of all the HSS/HSA wavelet 

filters which cannot be implemented by using solely WSA lifting steps. 

It therefore forms a solution to the problem of the implementation by 
lifting of pairs of filters (Hq,^) where the two filters have the same support and 
where the support is greater than 2. 
30 To this end, the present invention proposes a filtering method 

adapted to transform an input digital signal into one or more output digital 
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signals having even-indexed samples and odd-indexed samples, this method 
including at least one iteration which contains 

- an operation of modifying even-indexed samples by a function of 
weighted odd-indexed samples, 

5 - an operation of modifying odd-indexed samples by a function of 

weighted even-indexed samples, 

the weighted samples being obtained by at least one weighting 

operation, 

this method being remarkable in that at least one of the weighting operations is 
10 applied to the difference between two consecutive even-indexed samples. 

Thus, the present invention allows notably the implementation by 
lifting of any pairs of HSS/HSA wavelet filters, and the implementation by lifting 
of any pair of orthogonal filters and the implementation by lifting of any pair of 
any wavelet filters with equal support. 
15 The invention also makes it possible, by virtue of the technique of 

filtering by lifting, to reduce the number of filtering operations (multiplications 
and additions) compared with a conventional convolution filtering and also 
guarantees the reversible character of the wavelet transformations, that is to 
say a perfect reconstruction of the signal, in the case of signals whose sample 
20 values are integer. 

According to a particular characteristic, the operation of modifying 
odd-indexed samples is performed following the operation of modifying even- 
indexed samples. 

In a particular embodiment, the aforementioned iteration consists 

25 notably of: 

- weighting, by means of a first weighting coefficient, at least one 
odd-indexed sample adjacent to an even sample currently being modified, so as 
to obtain a weighted odd-indexed sample, 

- modifying at least one even-indexed sample using at least one 
30 weighted odd-indexed sample, 
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- weighting, by means of a second weighting coefficient, even- 
indexed samples adjacent to an odd sample currently being modified, so as to 
obtain weighted even-indexed samples, and 

- modifying at least one odd-indexed sample using at least one 
5 weighted even-indexed sample. 

"Adjacent" sample means a sample whose rank is consecutive with 
the sample in question. 

According to a particular characteristic, the second weighting 
coefficient is a function of the first weighting coefficient. 
10 The second weighting coefficient can for example depend on the first 

weighting coefficient as follows: 

f \ ^ 



Poj =r V 



i 

i=0 j 



where ao.i designates the first weighting coefficient, p 0l j designates the second 
weighting coefficient, i and j are integers and mj is a value defined by the 

15 recurrence m 0 =(-1) L ° and mj = -rrij_i, L 0 being a predetermined integer. 

According to a particular characteristic, at each iteration, the odd- 
indexed sample adjacent to the even sample currently being modified is 
alternately the sample of rank immediately lower or immediately higher. 

According to a particular characteristic, the filtering method includes, 
20 at the end of the aforementioned iteration, an additional filtering step including 
an operation of weighting by means of a third weighting coefficient. 

According to a particular characteristic, the third weighting coefficient 
is a function of the weighting coefficient used at the previous step, as follows: 

y = -i/(2Po lLo -i) 

25 where y designates the third weighting coefficient, L 0 is a predetermined 
parameter and Po,Lq-1 designates the weighting coefficient used at the 
previous step. 

In the application of the invention to the filtering of digital images, the 
input digital signal represents an image. 
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According to a second aspect, the present invention also proposes a 
filtering method adapted to transform one or more input digital signals into an 
output digital signal, the input signals including even-indexed samples and odd- 
indexed samples, the method including at least one iteration which contains 
5 - an operation of modifying odd-indexed samples by means of a 

function of weighted even-indexed samples, 

- an operation of modifying even-indexed samples by means of a 
function of weighted odd-indexed samples, 

the weighted samples being obtained by means of at least one 
10 weighting operation, 

this method being remarkable in that at least one of the weighting operations is 
applied to the difference between two consecutive even-indexed samples. 

According to a particular characteristic, the operation of modifying 
even-indexed samples is performed following the operation of modifying odd- 
1 5 indexed samples. 

In a particular embodiment, the aforementioned iteration consists 

notably of: 

- weighting, by means of a fourth weighting coefficient, even-indexed 
samples adjacent to an odd sample currently being modified, so as to obtain 

20 weighted even-indexed samples, 

- modifying at least one odd-indexed sample using at least one 
weighted even-indexed sample, 

- weighting, by means of a fifth weighting coefficient, at least one 
odd-indexed sample adjacent to an even sample currently being modified, so as 

25 to obtain a weighted odd-indexed sample, and 

- modifying at least one even-indexed sample using at least one 
weighted odd-indexed sample. 

According to a particular characteristic, the fourth weighting 
coefficient is a function of the fifth weighting coefficient. 
30 The fourth weighting coefficient can for example depend on the fifth 

weighting coefficient as follows: 
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Po.j =mj/ 



1 - 2 Z a o,i 
V i=o 



where a 0 j designates the fifth weighting coefficient, p 0j designates the fourth 
weighting coefficient, i and j are integer and m, is a value defined by the 
recurrence m 0 = (-1) L ° and nrij = -mj.-i, L 0 being a predetermined integer. 

According to the second aspect of the invention, according to a 
particular characteristic, at each iteration, the odd-indexed sample adjacent to 
the even sample currently being modified is alternately the sample of rank 
immediately lower or immediately higher. 

According to the second aspect of the invention, according to a 
particular characteristic, the filtering method includes, at the end of the 
aforementioned iteration, an additional filtering step including an operation of 
weighting by means of a sixth weighting coefficient. 

According to a particular characteristic, the sixth weighting coefficient 
is a function of the weighting coefficient used at the previous step, as follows: 

y=-i/(2Po,L 0 -i) 

where y designates the sixth weighting coefficient, l_ 0 is a predetermined 
parameter and Po,L 0 -i designates the weighting coefficient used at the 
previous step. 

In the application of the invention to the filtering of digital images, the 
output digital signal represents an image. 

According to a particular characteristic of the first and second 
aspects of the invention, the modification operations consist of applying an 
approximation function. 

This function may be: 

- the identity function, or 

- a function of a real variable which supplies the integer closest to 
the variable, or 

- a function of a real variable which supplies the first integer below 
the variable, or 
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- a function of a real variable which supplies the first integer above 
the variable, or 

- a function of a variable decomposed into subvariables whose sum 
is equal to the variable, which supplies a sum of approximate values of the 
subvariables, each of the approximate values of the subvariables being either a 
function of a real variable which supplies the integer closest to the variable, or a 
function of a real variable which supplies the first integer below the variable, or 
a function of a real variable which supplies the first integer above the variable. 

The invention also relates to a signal processing device which has 
means adapted to implement a filtering method as above. 

For the same purpose as indicated above, the present invention also 
proposes, according to a third aspect, a digital filtering device adapted to 
transform an input digital signal into one or more output digital signals 
containing even-indexed samples and odd-indexed samples, this filter having 

- at least one weighting module, 

- a module for modifying even-indexed samples by means of a 
function of weighted odd-indexed samples, 

- a module for modifying odd-indexed samples by means of a 
function of weighted even-indexed samples, 

these weighted samples being supplied by the weighting module, the 
modification modules functioning iteratively, so as to modify even-indexed 
samples at least once and then odd-indexed samples at least once, 
this filtering device being remarkable in that at least one of the aforementioned 
weighting modules receives as an input the difference between two consecutive 
even-indexed samples. 

For the same purpose as indicated above, the present invention also 
proposes, according to a fourth aspect, a digital filtering device adapted to 
transform one or more input digital signals into an output digital signal, the input 
signals containing even-indexed samples and odd-indexed samples, this 
filtering device having 

- at least one weighting module, 
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- a module for modifying odd-indexed samples by means of a 
function of weighted even-indexed samples, 

- a module for modifying even-indexed samples by means of a 
function of weighted odd-indexed samples, 

5 the weighted samples being supplied by the weighting module, the 

modification module functioning iteratively, so as to modify odd-indexed 
samples at least once and then even-indexed samples at least once, 
this filtering device being remarkable in that at least one of the aforementioned 
weighting modules receives as an input the difference between two consecutive 
10 even-indexed samples. 

The particular characteristics and the advantages of the filtering 
device being similar to those of the filtering method according to the invention, 
they are not stated here. 

Still with the same purpose, the invention also proposes a signal 
15 processing device including a filtering device as previously defined, or the 
means of implementing the corresponding method. 

The invention also concerns a signal processing device including at 
least two filtering devices as previously defined, the output signal of one of the 
filtering devices being the input signal of the other filtering device. 
20 The invention also concerns a digital apparatus including the signal 

processing device. 

The invention also concerns a digital photographic apparatus 
including the signal processing device. 

Still for the same purpose, the present invention also proposes an 
25 encoding method including steps adapted to implement a filtering method as 
above. 

Still for the same purpose, the present invention also proposes an 
encoding device including at least one filtering device as above, or the means of 
implementing the corresponding method. 
30 Still for the same purpose, the present invention also proposes a 

digital signal compression method, including steps adapted to implement a 
filtering method as above. 
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Still for the same purpose, the present invention also proposes a 
digital signal compression device, including at least one filtering device as 
above, or the means of implementing the corresponding method. 

An information storage means which can be read by a computer or 
5 by a microprocessor, integrated or not into the device, possibly removable, 
stores a program implementing a filtering method as above. 

The present invention also relates to a computer program product 
containing sequences of instructions for implementing a filtering method as 
above. 

10 As the particular characteristics and the advantages of the signal 

processing devices, of the digital apparatus, of the digital photographic 
apparatus, of the encoding and compression devices and methods, of the 
storage means and of the computer program product are similar to those of the 
filtering method according to the invention, they are not stated here. 

15 Other aspects and advantages of the invention will emerge from a 

reading of the following detailed description of a particular embodiment, given 
by way of non-limitative example. The description refers to the drawings which 
accompany it, in which: 

- Figure 1 depicts schematically a data processing device according 
20 to the invention, 

- Figure 2 depicts schematically another data processing device 
according to the invention, 

- Figure 3 depicts an embodiment of the data processing device 
according to the invention, 

25 - Figure 4 depicts an embodiment of a transformation circuit 

included in the data processing device of Figure 1 , 

- Figure 5 is a flow diagram illustrating the principal steps of a 
filtering method according to the present invention, in a first particular 
embodiment; 

30 - Figure 6 is a flow diagram illustrating the principal steps of a 

filtering method according to the present invention, in a second particular 
embodiment; 
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- Figure 7 is a flow diagram illustrating a method of carrying out the 
determination of the lifting parameters from the coefficients of the wavelet filters; 

- Figures 8A, 8B and 8C are flow diagrams illustrating respectively 
the calculation procedures denoted UNLIFTJHO, UNLIFTJH1 and UNLIFT_EQ 

5 acting in the flow diagram of Figure 7; 

- Figure 9 depicts a first embodiment of an elementary 
transformation unit according to the present invention; and 

- Figure 10 depicts a second embodiment of an elementary 
transformation unit according to the present invention. 

10 According to a chosen embodiment depicted in Figure 1, a data 

processing device according to the invention is a data encoding device 2 which 
has an input 24 to which a source 1 of non-encoded data is connected. The 
data processing device can be integrated into a digital apparatus, such as a 
computer, a digital photographic apparatus or a facsimile machine, for example. 

15 The source 1 has for example a memory means, such as a random 

access memory, hard disk, diskette or compact disc, for storing non-encoded 
data, this memory means being associated with an appropriate reading means 
for reading the data therein. A means for recording the data in the memory 
means can also be provided. The source 1 can also be integrated into the 

20 digital apparatus. 

It will be considered more particularly hereinafter that the data to be 
encoded are a series of digital samples representing an image IM. 

The source 1 supplies a digital image signal SI at the input of the 
encoding circuit 2. The image signal SI is a series of digital words, for example 

25 bytes. Each byte value represents a pixel of the image IM, here with 256 levels 
of grey, or black and white image. The image can be a multispectral image, for 
example a colour image having components in three frequency bands, of the 
red-green-blue or luminance and chrominance type. Each band is then 
processed in a similar manner to the monospectral image. 

30 Means 3 using encoded data are connected at the output 25 of the 

encoding device 2. 
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The user means 3 include for example means of storing encoded 
data, and/or means of transmitting encoded data. 

The encoding device 2 has conventionally, as from the input 24, a 
transformation circuit 21, to which the present invention more particularly 
5 relates, and an example embodiment of which will be detailed hereinafter. The 
transformations envisaged here are decompositions of the data signal into 
frequency sub-bands, so as to effect an analysis of the signal. 

The transformation circuit 21 is connected to a quantisation circuit 
22. The quantisation circuit implements a quantisation known per se f for 
10 example a scalar quantisation, or a vector quantisation, of the coefficients, or 
groups of coefficients, of the frequency sub-bands supplied by the circuit 21. 

The circuit 22 is connected to an entropic encoding circuit 23, which 
effects an entropic encoding, for example a Huffman encoding, or an arithmetic 
encoding, of the data quantised by the circuit 22. 
15 Figure 2 depicts another data processing device according to the 

invention, in the form of a device 5 for decoding data encoded by the device 2. 

Means 4 using encoded data are connected at the input 54 of the 
decoding device 5. The means 4 include for example means of storing encoded 
data, and/or means for receiving encoded data which are adapted to receive the 
20 encoded data transmitted by the transmission means 3. 

Means 6 using decoded data are connected at the output 55 of the 
decoding device 5. The user means 6 are for example image display means, or 
sound reproduction means, according to the nature of the data being 
processed. 

25 The decoding device 5 overall performs operations which are the 

reverse of those of the encoding device 2. The device 5 has an entropic 
decoding circuit 51, which performs an entropic decoding corresponding to the 
encoding of the circuit 23. The circuit 51 is connected to a dequantisation circuit 
52, corresponding to the quantisation circuit 22. 

30 The circuit 52 is connected to a reverse transformation circuit 53, 

corresponding to the transformation circuit 21. The present invention relates 
more particularly to the reverse transformation circuit 53. An example 
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embodiment will be detailed hereinafter. The transformations envisaged here 
effect a synthesis of the digital signal, from frequency sub-bands. 

The encoding device and the decoding device can be integrated into 
the same digital apparatus, for example a digital camera. In this case, the data 
5 processing device effects the encoding and the decoding of the data. 

With reference to Figure 3, an example of a device 10 implementing 
the invention is described. This device is adapted to transform a digital signal 
and, according to the examples developed hereinafter, to analyse it or to 
synthesise it or to analyse and then synthesise it. Thus Figure 3 illustrates an 
1 0 encoding or decoding device. 

The device 10 is here a microcomputer having a communication bus 
101 to which there are connected: 

- a central processing unit 100, 

- a read-only memory 102, 

15 - a random access memory 103, 

- a screen 104, 

- a keyboard 114, 

- a hard disk 108, 

- a disk drive 1 09 adapted to receive a diskette 1 1 0, 

20 - an interface 112 for communication with a communication network 

113, 

- an input/output card 106 connected to a microphone 111. 

The random access memory 103 stores data, variables and 
intermediate processing results, in memory registers bearing, in the description, 
25 the same names as the data whose values they store. The random access 
memory 103 contains notably: 

- a register "/c", containing the current value of the variable k defined 

below, 

- a register 7\ containing the current value of the variable j defined 

30 below, 

- a register "y/\ in which there are stored the values of the samples 
yi of the transformed image as they are calculated, 
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- a register "/3 k J\ in which there are stored the values of the 
weighting coefficients p k j defined below, 

- a register "m/\ in which there are stored the values of the 
parameters rrij defined below, 

5 - a register in which there are stored the values of the 

coefficients y defined below. 

In the case where Figure 3 illustrates an encoding device, the read- 
only memory 102 is adapted to store notably: 

- the operating program of the central processing unit 100, in a 
1 0 register "Program'', 

- the maximum value kmax of the variable k, in a register u k ma x\ 

- the value of the parameters U, defined below, for all the values of 
the variable k between 0 and kmax, in a register "L k \ 

- the values of the weighting coefficients a k ,j, in a register "a k J\ 

15 It should be noted that, in the case where Figure 3 illustrates a 

decoding device, the register program is also in read-only memory but the 
registers "k ma ", "U" and "a k f are in random access memory. 

The central processing unit 100 is adapted to implement the flow 
diagrams in Figures 5 and 7. 

20 The hard disk 108 stores the programs implementing the invention, 

as well as the data to be encoded and the encoded data according to the 
invention. These programs can also be read on the diskette 1 10, or received via 
the communication network 113, or stored in read-only memory 102. 

In more general terms, the programs according to the present 

25 invention are stored in a storage means. This storage means can be read by a 
computer or by a microprocessor. This storage means is integrated or not into 
the device, and may be removable. For example, it may include a magnetic 
tape, a diskette or a CD-ROM (fixed-memory compact disc). 

The device 10 can receive data to be encoded coming from a 

30 peripheral device 107, such as a digital photographic apparatus, or a scanner, 
or any other means of acquiring or storing data. 
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The device 10 can also receive data to be encoded coming from a 
distant device, via the communication network 113, and transmit encoded data 
to a distant device, again via the communication network 113. 

The device 10 can also receive data to be encoded coming from the 
5 microphone 111. These data are then a sound signal. 

The screen 104 enables a user notably to display the data to be 
encoded, and serves, with the keyboard 1 14, as a user interface. 

With reference to Figure 4, the transformation circuit 21, or analysts 
circuit, is a dyadic decomposition circuit with two levels. The circuit 21 is, in this 
10 embodiment, a conventional set of filters, respectively associated with 
decimators by two, which filter the image signal in two directions, respectively 
vertical and horizontal, into sub-band signals of high and low spatial 
frequencies. The relationship between a high-pass filter and a low-pass filter is 
determined by the conditions of perfect reconstruction of the signal. Different 
15 examples of filters will be envisaged hereinafter. It should be noted that the 
vertical and horizontal decomposition filters are not necessarily identical, 
although in practice this is generally the case. The circuit 21 has here two 
successive analysis units B1 and B2 for decomposing an image IM into sub- 
band signals according to two resolution levels. 
20 In generaf terms, the resolution of a signal is the number of samples 

per unit length used for representing this signal. In the case of an image signal, 
the resolution of a sub-band signal is related to the number of samples per unit 
length used for representing this sub-band signal horizontally and vertically. The 
resolution depends on the number of decimations carried out, on the decimation 
25 factor and on the resolution of the initial image. 

The first analysis unit B1 receives a digital image signal and applies it 
to digital filters, respectively low-pass and high-pass 210 and 220, which filter 
the image signal in a first direction, for example horizontal in the case of an 
image signal. After passing through decimators by two D210 and D220, the 
30 resulting filtered signals are respectively applied to two low-pass filters 230 and 
250 and high-pass filters 240 and 260, which filter them in a second direction, 
for example vertical in the case of an image signal. Each resulting filtered signal 
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passes through a respective decimator by two D230, D240, D250 and D260. 
The unit B1 outputs four sub-bands signals LLi, LHi, HLi and HHi with the 
highest resolution RES-i in the decomposition. 

The sub-band signal LLi contains components, or coefficients, of low 
5 frequency, in both directions, of the image signal. The sub-band signal LHi 
contains the components of low frequency in a first direction and of high 
frequency in a second direction, of the image signal. The sub-band signal HLi 
contains the components of high frequency in the first direction and the 
components of low frequency in the second direction. Finally, the sub-band 

10 signal HHi contains the components of high frequency in both directions. 

Each sub-band signal is a set of real coefficients constructed from 
the original image, which contains information corresponding to an orientation 
which is respectively vertical, horizontal and diagonal of the contours of the 
image, in a given frequency band. Each sub-band signal can be assimilated to 

15 an image. 

The sub-band signal LLi is analysed by the analysis unit B2 similar to 
the previous one in order to supply four sub-band signals LL 2 , LH 2 , HL 2 and HH 2 
of resolution level RES 2 . 

Each of the sub-band signals of resolution RES 2 also corresponds to 
20 an orientation in the image. 

To a given analysis circuit 21 there corresponds, in a conventional 
manner, a synthesis circuit whose structure is derived from that of the analysis 
circuit. 

Figure 5 shows the succession of lifting steps implementing the 
25 forward wavelet transformation according to the present invention. The reverse 
transformation will be described next in relation to Figure 6. 

An input signal x n to be transformed is considered, having samples 
whose ranks are between two integers i 0 and h (excluding h), into an output 
signal y n having samples whose ranks are between io and ii (excluding h). 
30 An initialisation step 500 consists first of all of setting the value of the 

samples of the output signal y n to that of the samples of the input signal x n . 
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During step 500, two counters, whose current value will be denoted 
respectively k and j, are also initialised to the value 0. 

During step 500, the parameter m jT mentioned above in relation to 

Figure 3, is also initialised, to a value m 0 ^(-1) L °, where L 0 is another 

5 parameter, defined below. 

If a signal y n is considered, having samples whose ranks are 
between i 0 and i 1 (excluding ii), the signal y n must be extended beyond this 
range, sufficiently for all the samples of the output signal y n whose ranks are 
between io and i 1 (excluding h) to be able to be calculated. This will emerge 
10 from the description of an example, given below. 

Thus a signal extension procedure (not depicted in Figure 5) consists 
of calculating the values of the output signal y for values of i leaving a 
predetermined interval io ^ i < H : 

Yi-ypsEOJo^) 04) 
1 5 where PSE(i,i 0 ,ii)) = io + min[mod(i-io,2(ii-io)),2(ii-io)-mod(Ho,2(i 1 -i 0 ))] (1 5) 

where min(a,b) designates the function supplying the smaller of the values a 
and b and mod(a,b) designates the integer lying between 0 and b-1 (including 
b-1) equal to the value of a, to within a multiple of b. 

The example, in no way limitative, of an extension procedure given 
20 above applies to signals such that the rows io and U are even. Other extension 
procedures can be used in the case where at least one of the rows io and U is 
odd. Nevertheless the remainder of the description applies to all parity cases. 
Moreover, there are other extension procedures for other types of filter, such as 
for example the so-called circular extension, known to a person skilled in the 
25 art, for orthogonal filters. 

The following step 502 of the filtering method is the first lifting step. 
This is a high-pass lifting step. 

In a preferred embodiment, this lifting step is either of the LSC or of 

the RSC type, according to the value of the parameter m 0 = (-1) L ° , and is 
30 defined by: 

Vn,y 2n+ 1 = Y2n + 1 -R(y2n + 1-m 0 ) ( 16 ) 
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This lifting step corresponds to the HSA high-pass filter H 1 (z) = z-1 
or 1 -z, according to the value of mo. it corresponds to an increase by one unit 
of the support of the equivalent high-pass filter. 

At the end of the lifting step 502, a test 504 is carried out in order to 
5 determine whether or not the value of the counter j is equal to L 0 . 

As long as the result of the test 504 is negative, a step 506 is 
iterated, which consists essentially of a sequence of operations of modifying 
and weighting the even-indexed or odd-indexed samples of the signal y n . 

This is because, during the step 506, a succession of L 0 pairs of 
10 lifting steps are effected, corresponding to values of the counter j between 0 
and 0 and L 0 -1, where the first lifting step is a single coefficient low-pass lifting 
step: 

Vn,y 2n =y2n + R(ao i j-y2n+m J ) ( 17 ) 
and where the second lifting step is a HSA high-pass lifting step: 
15 Vn,y 2n+1 =y 2n+ i +R(Po,j-(y2n ~V2n+2)) (18) 

The non-zero parameter a 0 j used in equation (17) is a function of 

the coefficients of the pair of wavelet filters (Hq,^) and its calculation is 
detailed below. 

Given that the value of mj is either -1 or +1 as will be seen, the rank 
20 of the sample y2n+mj is either 2n-1 or 2n+1 , that is to say an integer following 

on from the rank of the sample y 2n currently being modified. It is said that the 
odd-indexed sample y2n+mj is adjacent to the sample y 2n . 

The single coefficient low-pass lifting step is either of the LSC or 
RSC type, according to the value of the parameter mj, which is defined 
25 recursively by: mj = -m^ . 

There is a systematic alternation of the following single coefficient 
low-pass lifting steps (that is to say corresponding to following values of j): for 
example, the second is of the LSC type if the first is of the RSC type or vice 
versa. The last single coefficient low-pass lifting step is always of the LSC type. 
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The low-pass lifting step of equation (17) guarantees that the support 
of the equivalent low-pass filter increases until it reaches the value of the 
support of the equivalent high-pass filter of the previous high-pass lifting step, 
i.e. the closest higher even value. 
5 In accordance with a preferred embodiment of the present invention, 

the high-pass lifting step of equation (18) is of the HSA type, unlike the prior art, 
which uses solely single coefficient lifting steps and WSA lifting steps for the 
implementation by lifting of HSS/HSA wavelet filters. 

The high-pass lifting step of equation (18) guarantees that the 
10 support of the equivalent high-pass filter increases by a value of two (that is to 
say until it reaches the closest higher even value) with respect to the support of 
the equivalent high-pass filter of the previous high-pass lifting step. 

The equivalent high-pass filter of the high-pass lifting step of 
equation (18) is of the HSA type if and only if the equivalent high-pass filter of 
15 the previous high-pass lifting step is HSA and the non-zero parameter p 0) j use d 

in equation (18) is a function of the parameter mj and of the parameters a 0 i 

originating from the previous single coefficient low-pass lifting steps, this 
function being for example defined as follows: 

f i ^ 



Poj = m j ' 



!=0 j 



(19) 



20 Unlike the low-pass lifting step of equation (17), which lifts the low- 

pass coefficients from the high-pass coefficients alternately to left and right, the 
high-pass lifting step of equation (18) lifts the high-pass coefficients from both 
the left-side and right-side low-pass coefficients. 

When the result of test 504 is positive, that is to say j = Lo, the 

25 alternation of SC low-pass and HSA high-pass lifting steps ends and is followed 
by an RSC low-pass lifting step 508, defined by: 

Vn, y 2n = y 2n + R(y-Y 2n+1 ) (20) 
In a preferred embodiment, the coefficient y is calculated according to 
the weighting coefficient used at the previous high-pass lifting step: 
30 Y = -1/(2Po b l 0 -i) - (21) 
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in order to guarantee that the corresponding low-pass filter is HSS. 

The low-pass lifting step 508 guarantees that the support of the 
equivalent low-pass filter increases until it reaches the value of the support of 
the equivalent high-pass filter of the previous high-pass lifting step, that is to 
5 say, in the preferred embodiment described here, the closest higher even 
integer. 

In this preferred embodiment, there is obtained, at the end of step 
508, a corresponding pair of low-pass and high-pass filters which have equal 
even supports and which are respectively HSS and HSA. 
10 Next an iteration is carried out of steps of WSA low-pass (step 512) 

and WSA high-pass (step 518) lifting, that is to say k max - WSA lifting steps in 
total, the value of the counter k increasing from 1 to k max - . 

For this purpose, following step 508, the value of a counter k is 
incremented by one unit at step 510. The value of k is incremented once again 
15 after each passage through step 512 (see step 514 in Figure 5) and after each 
passage through step 518 (see step 520 in Figure 5). 

Steps 512, 514, 518 and 520 are reiterated until the value of k is 
equal to k ma x (see tests 516 and 522 on the value of k, respectively carried out 
after the incrementation steps 514 and 520). 
20 The WSA low-pass lifting steps 512 are defined by: 

-2k-r 



Vn,y 2n =Y2n+R 



X a 2k-1,j-(y2n+1-2j -Y2n+1+2j) 
J=0 



(22) 



and the WSA high-pass lifting steps 518 are defined by 



Vn,y 2n+ 1 =Y2n + 1 



X!^2k t j-(y2n-2j ~Y2n+2j) 



(23) 



These steps are known per se. Their effect is an increase in the size 
25 of the corresponding filter (low-pass for step 512 and high-pass for step 518) by 
a multiple of four. 

In order to produce wavelet filters with a unity direct current gain and 
a Nyquist gain of 2, that is to say with a so-called (1,2) normalisation, no 
specific normalisation step of the type y» = G.yj, where G is a real value, is 
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necessary. On the other hand, if a different normalisation is required, such as a 
(Gi, 2G 2 ) normalisation with a DC gain of Gi and a Nyquist gain of 2G 2) a final 
normalisation step is necessary: y 2n =Gi-V2n anc ' Y2n+1 = G 2 .y2n+i* 

The approximation function R is not necessarily the same for all the 
5 lifting steps for which it is used. 

The parameters to be chosen for the lifting implementation of the 
HSS/HSA filters are Lk f a^j,k max and the function R at each lifting step, it is 

thus possible to establish a one-to-one mapping between these parameters and 
the coefficients and lengths of any pair of perfect reconstruction HSS/HSA 
10 wavelet filters. 

Figure 6 shows the succession of lifting steps implementing the 
reverse wavelet transformation in accordance with the present invention. 

An input signal y n to be transformed is considered, having samples 
whose ranks lie between two integers i 0 and ii (not including U) and an output 
15 signal x n having samples whose ranks lie between i 0 and it (not including h). 

An initialisation step 600 consists first of all of setting the value of the 
samples of the output signal x n to that of the samples of the input signal y n . 

All the steps of the reverse lifting implementation are the operations 
which are the reverse of the lifting steps performed in forward transformation. 
20 Moreover, the steps of the reverse lifting implementation are performed in the 
opposite order to the steps of the forward lifting implementation. 

For an analysis processing followed by a synthesis processing of the 
same signals, all the parameters of the reverse transformation are identical to 
those chosen for the forward transformation. 
25 As shown by Figure 6, during the initialisation step 600, a counter, 

whose current value is denoted k, is also initialised to the value k max - 1 . 

If a signal x n is considered, having samples whose ranks lie between 
i 0 and h (not including ii), the signal x n must be extended beyond this range, 
sufficiently for all the samples of the output signal x n whose ranks are between 
30 i 0 and h (not including h) being able to be calculated. This will emerge from the 
description of an example, given below. 
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Thus a signal extension procedure (not depicted in Figure 6) consists 
of calculating the values of the output signal Xj for values of i leaving a 
predetermined interval i 0 < i < ii , 

- as defined by equation (24) below for values of i such that mod(i- 
5 i 0 ,2(ii-i 0 )) < Mo : 

x i = x io+mod(i-i 0i 2(i 1 -io)) ( 24 ) 
where mod(a,b) designates the integer lying between 0 and b-1 (including b-1) 
equal to the value of a, to within a multiple of b, 

- as defined by equation (25) below for even values of i such that 
1 0 mod(H 0l 2(Mo)) ^ Mo : 

x i = x i 0 +mod(2(i 1 -i 0 -1)-(i-i 0 ),2<H -\ 0 ))) (25) 

- and as defined by equation (26) below for odd values of i such that 
mod(i-i 0 ,2(Mo)) £ Mo : 

x i = x i 0 +mod(2(i 1 -i 0 )-(M 0 -i 0 ))) < 26 > 
1 5 There are other extension procedures for other types of filter, such as 

for example the so-called circular extension, known to a person skilled in the 
art, for orthogonal filters. 

In order to implement wavelet filters with a unity DC gain and a 
Nyquist gain of 2, that is to say with a so-called (1,2) normalisation, no specific 
20 normalisation step of the type Xj = G.x,, where G is a real value, is necessary. 
On the other hand, if a different normalisation is required, such as a (Gi, 2G 2 ) 
normalisation with a DC gain of Gi and a Nyquist gain of 2G2, a reverse 
normalisation step is necessary: x 2n =(1/G 1 ).x 2n and x 2n+ i =(1/G2)- X 2n+1- 

As shown in Figure 6, an iteration of steps of WSA low-pass lifting 
25 (step 602) and WSA high-pass lifting (step 608) is carried out, that is to say 
k max - 1 WSA lifting steps in total, the value of the counter k decreasing from 
kmax -1 to1. 

For this purpose, the value of the counter k is decremented by one 
unit after each passage through step 602 (see step 604 in Figure 6) and after 
30 each passage through step 608 (see step 610 in Figure 6). 
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Steps 602, 604, 608 and 610 are reiterated until the value of k is 
equal to 1 (see tests 606 and 612 on the value of k, respectively performed after 
the decrementation steps 604 and 610). 

The WSA low-pass lifting steps 602 are defined by: 



Vn,x 2n =x 2n -R 



L 2k-1 - 1 

Z! a 2k-1,j-( x 2n+1-2j -X 2n+ i + 2j) 



(27) 



and the WSA high-pass lifting steps 608 are defined by: 

fL 2k -1 ^ 



Vn,x 2n+1 =x 2n+ i-R 



Z a 2k,j-( x 2n-2j ~ x 2n+2j) 



(28) 



These steps are known per se. 

When the result of test 606 or 612 is positive, that is to say k = 1 , the 
10 alternation of the steps of WSA low-pass lifting and WSA high-pass lifting ends 

and is followed by a step of RSC low-pass lifting 614, defined by: 

Vn,x 2n = x 2n -R(y.x 2n+1 ) (29) 
In a preferred embodiment, the coefficient y is calculated as a 

function of the weighting coefficient used at the following high-pass lifting step: 

15 Y = -1/(2p 0l Lo-l) (30) 

During step 614, the parameter mj is also initialised to the value 1 
and a counter] to the value I_o-1 . 

At the end of the low-pass lifting step 614, a test 616 is performed in 
order to determine whether or not the value of the counter j is equal to 0. 
20 As long as the result of test 616 is negative, a step 618 is iterated, 

which consists essentially of a sequence of operations of modifying and 
weighting the even-indexed or odd-indexed samples of the signal x n . 

This is because, during step 618, a succession of I_o pairs of lifting 
steps are performed, corresponding to values of a counter j lying between Lo-1 
25 and 0, where the first lifting step is a HSA high-pass lifting step: 

Vn ^ 2n+1 - x 2n+1 -R(p 0| j-(X2n -x 2n+2 )) (31) 
and where the second lifting step is a single coefficient low-pass lifting step: 

Vn,x 2n - x 2n -R(a 0 j.x 2n+mj ) (32) 
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When the result of test 616 is positive, that is to say j = 0, the last 
lifting step 620 is performed, which is the reverse of the first high-pass lifting 
step performed in the forward transformation. In a preferred embodiment, this 
lifting step is either of the LSC or of the RSC type, according to the value of the 

5 parameter m 0 = (-1) L ° , and is defined by: 

Vn,x 2n+1 =x 2n+1 +R(x 2n+1 _ mo ) (33) 

Figure 7 illustrates a procedure for determining the lifting parameters 
for any given pair of HSS/HSA wavelet filters, defined by equations (4), (6) and 
(7) appearing in the introduction. This procedure also applies to orthogonal 
10 filters as well as to any filters, without any particular symmetry, with equal 
support. The "UNLIFT" procedures appearing below are detailed subsequently. 

As shown by Figure 7, a first step 700 consists of initialising an 
integer variable k to 0 and a variable no to -1 and storing, in a variable Ij, i e 
{0,1}, the length, that is to say the number of coefficients, of the intermediate 
15 low-pass and high-pass filters H 0 and Hi, respectively (that is to say the number 
of coefficients of the so-cafled equivalent filters). 

Then a step 702 consists of testing whether l 0 < h. 
If the result of test 702 is positive, there is performed a step 704 of 
determining the coefficients of the WSA high-pass lifting step 518 of Figure 5 by 
20 applying a procedure UNLIFT_H1 (see Figure 8B), at the end of which there are 
obtained the number of lifting coefficients U, the L_k coefficients for 0 < i < U 
and the new intermediate high-pass filter Hl 

At step 704, the value of k is also incremented by one unit, the length 
li of the intermediate filter is updated and the value 1 is stored in the variable no. 
25 If the result of test 702 is negative, a step 706 is performed, 

consisting of testing whether l 0 > h. 

If the result of test 706 is positive, a step 708 of determining the 
coefficients of the WSA low-pass lifting step 512 of Figure 5 is performed by 
applying a procedure UNLIFT J-I0 (see Figure 8A), at the end of which there are 
30 obtained the number of lifting coefficients L kj the Lr coefficients ock.i for 0 < i < L k 
and the new intermediate low-pass filter H 0 . 
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At step 708, the value of k is also incremented by one unit, the length 
li of the intermediate filter is updated and the value 0 is stored in the variable n 0 . 

If the result of test 706 is negative, a step 710 is performed, 
consisting of testing whether the variable n 0 is zero. 
5 If the result of test 710 is positive, a step 712 is performed, consisting 

of storing the value 0 in the variable L k and incrementing the value of k by one 
unit 

At the end of step 712, or if the result of test 710 is negative, a step 
714 of determining the coefficients of the single coefficient low-pass lifting step 
10 included in step 506 in Figure 5 is performed by applying a procedure 
UNLIFTJEQ (see Figure 8C), at the end of which the number of lifting 
coefficients L 0 and the L 0 coefficients a 0l i for 0 < i < L 0 are obtained. 

At step 714, the value of k is also incremented by one unit. 
At the end of step 714, a rearrangement step 716 is performed, 
15 consisting of updating the lifting coefficients obtained by reordering them, as 
follows: 

(Lo,Li,...,Lk_i) = (l_k-i,Lk-2,.'ML-o) 
(ao,i,ai i j,...,ak-ij) = (ak-i^ak^b.^ccoj) 
During step 716, the value k max is also stored in the variable k. 
20 In summary: 

- when lo > li, the lifting coefficients are obtained by the procedure 
UNLIFT_H0(H 0 ,Hi) f 

- when l 0 < U, the lifting coefficients are obtained by the procedure 
UNLIFT_H1(Ho,Hi) and 

25 - when l 0 = h, the lifting coefficients are obtained by the procedure 

UNLIFT_EQ(H 0f Hi). 

These three procedures are now described in more detail. 
As shown by Figure 8A, which illustrates the procedure UNLIFT_H0, 
an initialisation step 800A consists of storing in the variable l 0 half of the length 
30 of Ho and storing in the variable U half of the length of Hi. During step 800A, a 
variable i and the variable Lk are also initialised to the value 0. 
Then a step 802A consists of testing whether l 0 > 



29 



If the result of test 802A is negative, the procedure UNL1FT_HG is 
terminated. 

If the result of test 802A is positive, a step 804A consists of 
calculating and H 0 (z), as follows: 

a k,i = ~ H 0,l 0 /H 1,H 

H 0 (z) = H 0 (z) - a kJ .(z-' 0 +1 1 - z ! ° -'1 ^ (z) 

At the end of step 804A, the value of l 0 is decremented by two units 
and the value of i and the value of Lk are incremented by one unit. 
Then this step 804A is reiterated as long as l 0 > h. 
As shown by Figure 8B, which illustrates the procedure UNLIFT_H1 f 
10 an initialisation step 800B identical to step 800A of Figure 8A is first of all 
performed. 

Then a step 802B consists of testing whether U > l 0 . 
If the result of test 802B is negative, the procedure UNLIFTJ-M is 
terminated. 

15 If the result of test 802B is positive, a step 804B consists of 

calculating a k ,i and H<i(z), as follows: 

a k,i =-Hl i l 1 /H 0J 0 

^ (z) = Hi (z) - a kJ .(Z~ ! 1 +l o - Z ! 1 - { ° ).H 0 (z) 

At the end of step 804B, the value of h is decremented by two units 
and the value of i and the value of Lk are incremented by one unit. 
20 Then this step 804B is reiterated as long as l 1 > l 0 . 

As shown in Figure 8C, which illustrates the procedure UNLIFT_EQ, 
an initialisation step 800C consists of storing in the variable l 0 half of the length 
of H 0 , and storing in the variable h the value of l 0 . During step 800C, the 
variable L 0 is also initialised to the value I 0 -1 and a variable i to the value L 0 -1 . 
25 The following step 802C consists of calculating Ho(z) and the 

parameter y, as follows: 

Y = H o,l 0 

H 0 (z) = z.[H 0 (z)-y.H 1 (z)] 
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At the end of step 802C, the value of i 0 is decremented by one unit 
and a variable q is initialised to the value 0. 

Then a step 804C consists of testing whether lo > 0. 

If the result of test 804C is negative, the procedure UNLIFTJEQ is 
terminated. 

If the result of test 804C is positive, a step 806C consists of testing 
whether q > 0. 

If the result of test 806C is positive, the following step 808C consists 
of calculating p 0 j. H-i(z), oto,i and Ho(z), as follows: 

Poj = H 1 f -li+1 /H 0 I -lo+'I 

H 1 (z) = z- 1 -IH 1 (z)-Po i ,.(z- 1 -z)Wo(2)] 

a 0,i =~ H OHo+1 /H Vli+1 
H 0 (z) = z- 1 .[Ho(z)-a 0 i.H 1 (z)] 



During step 808C, after calculating H«i(z), the value of h is 
decremented by one unit and, after calculating H 0 (z), the value of l 0 and the 
value of i are decremented by one unit and the value 1 is stored in the variable 

q- 

15 Then this step 808C is reiterated as long as l 0 > 0. 

If the result of test 806C is negative, the following step 810C consists 
of calculating p 0 j, H-i(z), a 0 ,i and Ho(z), as follows: 

Po,i =-*" , 1 ( i 1 /H 0J 0 

H 1 (z) = z.[H 1 (z)-p 0j i-(2" 1 "Z).Ho(z)] 

a 0,i = H 0,I 0 /H Ui 

H 0 (z) = z.[H 0 (z)-a oj .H 1 (z)] 

During step 81 0C, after calculation of H^z), the value of h is 
20 decremented by one unit and, after calculation of Ho(z), the value of lo and the 
value of i are decremented by one unit and the value 0 is stored in the variable 

q. 

Then this step 81 0C is reiterated as long as l 0 > 0. 
Figures 9 and 10 illustrate in schematic form an elementary 
25 transformation unit included in a filtering device in accordance with the present 



31 



invention, in two particular embodiments. Figure 9 illustrates more particularly a 
forward lifting transformation unit and Figure 10 illustrates more particularly a 
reverse lifting transformation unit. 

The filtering device 900 shown schematically in Figure 9 is equivalent 
5 to that of the unit BE illustrated in Figure 4. 

It has an input E 20 to which the signal to be transformed is applied. 
The signal to be transformed includes a series of samples {Xj}, where i is a 
sample rank index. 

The input E20 is connected to a first decimator by two D 2 o, which 
10 delivers the even-indexed samples {x 2 J. 

The input E 20 is also connected to an advance AV21 followed by a 
second decimator D21, which delivers the odd-indexed samples {x 2i+ i}. 

The first decimator D 20 is connected to a first filter A 0 , which filters the 
even-indexed samples by applying to them a weighting coefficient p 0 j and an 
15 approximation function as described in equation (18) and supplies them to an 
adder AD 20 . 

The assembly formed by the filter A 0 and the adder AD 2 o is a high- 
pass lifting module, denoted FHP in the drawing. 

The output of the second decimator D21 is also connected to the 
20 adder AD 20 . 

Consequently the adder AD20 delivers intermediate samples y 2 t+i 
calculated according to equation (18). 

The output of the adder AD 2 o is also connected to a second filter A1 
which filters the odd-indexed samples by applying to them a weighting 
25 coefficient a 0) j and an approximation function as described in equation (17) and 
supplies them to an adder AD21. 

The assembly formed by the filter A1 and the adder AD 2 i is a low- 
pass lifting module, denoted FLP in the drawing. 

The output of the first decimator D 2 o is also connected to the adder 

30 AD21. 
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The output of the adder AD 21 is the first output S 2 o of the 
transformation unit, which delivers a digital signal containing low-frequency 
samples y 2i . 

The adder AD 2 o is connected to a second output S21 of the 
5 elementary unit, which delivers a digital signal containing high-frequency 
samples y 2 i+i. 

The filtering device 1000 shown schematically in Figure 10 effects 

the transformation which is the reverse of that effected by the filtering device 

900 shown schematically in Figure 9. 
10 This transformation unit has a first input E30 to which a first signal to 

be transformed is applied and a second input E 31 to which a second signal to be 

transformed is applied. 

The signals to be transformed contain here the samples obtained 

after analysis filtering of a digital signal by the analysis unit of Figure 9. These 
15 samples have possibly been modified by another processing between analysis 

and synthesis. 

More precisely, the first signal to be transformed contains the low- 
frequency samples {y 2 i} and the second signal to be transformed contains the 
high-frequency samples {y 2 j+i}. 
20 The synthesis unit has a structure similar to that of the analysis unit, 

and can be derived from it simply. In particular, the synthesis unit uses the 
same filters A 0 and as the analysis unit. 

The input E 30 is connected to a subtracter SO 30 . 

The input E 31 is connected to the filter Ai f itself connected to the 
25 subtracter S0 3 q. 

The assembly formed by the filter A1 and the subtracter SO30 is a 
low-pass lifting module, denoted ILP in the drawing. 

The output of the subtracter SO30 delivers a digital signal containing 
intermediate samples y 2i calculated according to equation (32). 
30 The output of the subtracter SO30 is connected to the filter Ao, itself 

connected to a subtracter SO31, to which the input E31 is also connected. 
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The assembly formed by the filter Ao and the subtracter SO31 is a 
high-pass lifting module, denoted IHP in the drawing. 

The output of the subtracter S0 31 delivers a digital signal containing 
reconstructed odd-indexed samples x 2i+ i. 
5 The output of the subtracter SO 30 delivers a digital signal containing 

reconstructed even-indexed samples x 2 i. 

The output of the subtracter SO30 is connected to an interpolator by 
two IN30 and the output of the subtracter S0 31 is connected to an interpolator by 
two IN31 itself connected to a delay R31. The interpolator IN 30 and the delay R 31 
10 are connected to an adder AD 32 , which delivers as an output S 30 the signal 
containing the reconstructed samples x t . 

An example, in no way limitative, of forward and reverse 
transformation by lifting for the implementation of wavelet filters according to the 
invention is now given. 
1 5 The pair of filters used in this example are as follows: 

u , x OC.Z" 1 +1 + Z + C6.Z 2 , x 

Hr)(z) = (35) 

0V J 2.(cc + 1) 

u , x a.z" 1 + 1-z + a.z 2 

Hl (z) = 57 T\ (36) 

2.(a-1) 

The lifting parameters corresponding to this pair of filters are: 
k = 1 , L 0 = 1 and a 0 ,o = a. (37) 
20 An approximation function R identical for all the lifting steps, defined 

by R(x) = E(x+1/2) where E designates the integer part, is also chosen. 

The length of the signal is chosen so as to be equal to 6: io = 0 and ii 

= 6. 

With regard to the forward lifting equations, all the values of the 
25 output signal Y are initialised to the values of the input signal X. The signal Y is 
then extended by one sample at each limit: y_i = y 0 and y 6 = y$- 

The following four lifting operations are performed. First of all, the 
values of the samples (y.i, y 1f y 3 , y 5 ) are calculated from: 

y2n + i = y 2 n + i-E(y 2n+ 2+1/2) for n = -1 , 0, 1 , 2 (38) 
30 Next, the values of the samples (y 0f y 2t y4, ye) are calculated from: 
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where 



where 



Y2n = y2n+E(a 0l o.y 2 n-i+1/2) for n = 0, 1 , 2, 3 (39) 

oto.o = a (40) 
Then the values of the samples (yi, y3, y$) are calculated from: 

y2n + 1 = y2n + 1 + E[po,0-(y2n-y2n + 2) + 1/2)] for n = 0, 1 , 2 (41 ) 

p 0l o = -1/(1-2a) (42) 
Finally, the values of the samples (y 0 , y 2 , y4) are calculated from: 

y 2 n = y2n+E(y.y 2 n + i+1/2) for n = 0, 1 , 2 (43) 



1 0 where 



Y = -1/(2Po.o) (44) 
With regard to the reverse lifting equations, all the values of the 
output signal X are initialised to the values of the input signal Y. The signal X is 
then extended by two samples at each limit: x_ 2 = x 0 , x_i = -xi and x 6 = X4 and xy 
1 5 = -x 5 . 

The following four reverse lifting operations are performed. First of 
all, the values of the samples (x. 2 , x 0 , x 2 , X4, x 6 ) are calculated from: 

x 2n = x 2n -E(y.x 2n+1 +1/2) for n = -1 , 0, 1 , 2, 3 (45) 
Next, the values of the samples (x_i, xi, X3, X5) are calculated from: 
20 x 2n+1 = x 2n+ rE[p 0> o.(x 2 n-x 2 n +2 )+l/2)] for n = -1 , 0, 1 , 2 (46) 

Then the values of the samples (x 0 , x 2l X4, X6) are calculated from: 
x 2n = x 2n -E(a 0 ,o.x 2n -i+1/2) for n = 0, 1 , 2, 3 (47) 
Finally, the values of the samples (xi, X3, X5) are calculated from: 
x 2n+ i = x 2n+ i+E(x 2n+2 +1/2) for n = 0, 1 , 2 (48) 
25 The field of application of the present invention is not limited to 

HSS/HSA wavelet filters with equal supports, but extends much more broadly, 
on the one hand, to orthogonal wavelet filters and, on the other hand, to any 
wavelet filters with equal supports. 

With regard to the application of the invention to the implementation 
30 by lifting of orthogonal wavelet filters, the only difference compared with what 
was described previously lies in the value of the lifting parameters aoj, poj. y and 
in the relationship between these parameters, in order to guarantee that they 
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implement HSS/HSA filters or orthogonal filters equally well. It should be noted 
that the method of calculating the lifting parameters, described above in relation 
to Figures 7 and 8, can be used for generating all the lifting parameters a 0tU p 0 ,j, 
y for the HSS/HSA filters and for orthogonal filters. 
5 All the orthogonal filters have equal support and even length but are 

not symmetrical, except in one trivial case: the Haar filter, defined by H 0 (z) = 
(1+z)/2 and H 1 (z) = -1+z. 

The orthogonal filters of length 4 are characterised by a parameter 0 
representing one degree of freedom: 

-1 2 
Ar . cos6.(sinG + cosO).z + sin0.(sin0 + cos 6) + sine.fstnO - cos6).z + cosG.(cos0 - sin0).z 

10 h 0 (z) = ' - — (49) 

2 

H^z) = cose. (cose - sine).z - sine.(sine - cos 9) + sin0.(sine + cos6).z - cose.(stne + cos8).z (50) 

In more general terms, the orthogonal filters satisfy the equation: 

H 1 (z)-2z.H 0 (-z" 1 ) (51) 
In the non-limitative example given here, as before, an approximation 
15 function R identical for all the lifting steps is chosen and defined by R(x) = 
x+1/2. 

The length of the signal is chosen so as to be equal to 6 : i<> = 0 and ii 

= 6. 

The implementation by lifting of such a filter is carried out using the 
20 same steps as in the previous example (except with regard to the extension 
procedure, described below) except that the parameters are here chosen as 
follows: 

k = 1 , L 0 = 1 , a 0 o = CQSe , p 0 o = cos(26), y = $ine ~ CQSe (52) 
0lU sine + cose p ' 2(sine + cos9) v ; 

The lifting factorisation procedure given above is similar for 
25 orthogonal filters. The only difference lies in the relationship between the 
coefficients a 0 j, Po.i, y. 

As has been seen for HSS/HSA filters, the parameters ao,i are a 
function of the coefficients of the filters of the wavelet transformation to be 
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effected. The same applies to orthogonal filters. The method of calculating the 
lifting parameters described above is unchanged. 

As has been seen for the HSS/HSA filters, the parameters p 0 ,i are a 
function of the parameters ocoj and are such that they guarantee that the 
5 equivalent high-pass filter Ht(z) of each high-pass lifting step is HSA. The 
situation is similar for orthogonal filters: the parameters poj are a function of the 
parameters ocoj and are such that they guarantee that the equivalent high-pass 

filter Hi(z) of each high-pass lifting step is such that H 0 (z) = -z.H 1 (-z~ 1 )/2 and 
Hi(z) define an orthogonal wavelet transformation. The method of calculating 
10 the lifting parameters described above can also be used for calculating the 
parameters po,i- 

Finally, as has been seen for the HSS/HSA filters, the parameter y is 
a function of the parameters Po,l 0 -i ancl * s suc ^ that jt guarantees that the 

equivalent low-pass filter H 0 (z) of this low-pass lifting step is HSS. The situation 

15 is similar for the orthogonal filters: the parameter y is a function of the 
parameters ocoj and p 0 j and is such that it guarantees that the equivalent low- 
pass filter H 0 (z) of this low-pass lifting step is such that H 0 (z) = -z.H 1 (-z" 1 )/2. 
The method of calculating the lifting parameters described above can also be 
used for calculating the parameter y. 

20 For transformation by forward lifting, all the values of the output 

signal Y are initialised to the values of the input signal X. The signal Y is then 
extended by one sample at each limit: in the example supplied here, use is 
made of an extension procedure known to a person skilled in the art as circular 
extension: y.-( = y 5 and ye = yo* 

25 All the subsequent steps are identical to those of the previous 

example, using the values given in equation (52). 

For transformation by reverse lifting, all the values of the output 
signal X are initialised to the values of the input signal Y. The signal X is then 
extended by two samples at each limit: x_2 = X4, x_i = X5, X6 = xo and Xr = xi. 

30 All the subsequent steps are identical to those of the previous 

example, using the values given in equation (52): 
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With regard to the application of the invention to the implementation 
by lifting of any wavelet filters with equal support, the only differences compared 
with what was described above for the HSS/HSA filters and the orthogonal 
filters are as follows: 

- there is no dependency relationship between the parameters ocqj, 
Po,i> y ; it is therefore not necessary to calculate these values in the RAM; 

- the first high-pass lifting step in the forward implementation by 
lifting can be any, that is to say: 



Y2n+1 == y2n+1 " R 



n+2i 



(53) 



10 where §j is a weighting coefficient. 

It should be noted that, for example, if use is made, in the forward 
transformation, instead of equation (16), of the following equation: 

y 2 n + i-y 2 n + i-R( y2n+ 2 y2n+2 ) d6a) 

then, in the reverse transformation, the last lifting step 620 is exactly its reverse, 
1 5 that is to say, in the above case: 

x 2 n + l-x 2n+1+ R[ X2n+ 2 X2n+2 ] (34) 

All the lifting coefficients ao,\, po.i, Y necessary for the implementation 
by lifting of a pair of any wavelet filters with equal support can be obtained by 
means of the method of calculating the lifting parameters described above in 
20 relation to Figures 7 and 8. The lifting coefficients 5j are then the even-indexed 
coefficients of the filter Hi(z) obtained from the formula (Lk, ak,i, Hi, H 0 ) = 
UNLIFT_EQ(H 0 ,Hi) of step 714 illustrated in Figure 7. 



